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We investigate tiie non-Abelian Josepiison effect in spinor Bose-Einstein condensates with double 
optical traps. We propose, for the first time, a real physical system which contains non-Abelian 
Josephson effects. The collective modes of this weak coupling system have very different density 
and spin tunneling characters comparing to the Abelian case. We calculate the frequencies of the 
pseudo Goldstone modes in different phases between two traps respectively, which are a crucial 
feature of the non-Abelian Josephson effects. We also give an experimental protocol to observe this 
novel effect in future experiments. 
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Introduction. — The Josephson effect is a quantum tun- 
neling phenomena which occurs when a pair of supercon- 
ducting or superfluid systems are connected together by 
a weak link due to some kind of physical barrier. Be- 
ginning with Josephson's original paper in 1962 the 
Josephson effect has become a paradigm of the phase co- 
herence manifestation in a macroscopic quantum system. 
With the amazing experimental progress in cold atom, 
the Josephson junction has been reahzed for the trapped 
Bose-Einstein condensates (BEC) of Ql- However, 
most of extensive studies about this effect focus on the 
Abelian case so far, in terms of a junction of two systems 
with spontaneously broken Abehan symmetry |j, @ . 

Recently, F. Esposito et al. generalized the Abelian 
Josephson effect to the non-Abelian case in a field the- 
oretic language The non-Abelian Josephson effect 
emerges in a junction of two systems with spontaneously 
broken non-Abelian gauge symmetries, which are weakly 
coupled together. These two systems should have the 
same symmetry so that there is an initial doubling of 
symmetry. This doubled symmetry should be sponta- 
neously broken and leads to a doubling of the corre- 
sponding Goldstone modes. However, due to the coupling 
between the two systems, only the diagonal part of the 
symmetry group generators is preserved, which will give 
rise to one set of true gapless Goldstone modes and one 
set of pseudo Goldstone modes. These pseudo Goldstone 
modes is gapped and leads to the non-Abelian Josephson 
effect. The non-Abelian Josephson effect is an universal 
effect in physics and it exists in a lot of physical system 
with a non-Abelian symmetry. There are some other pos- 
sible scenarios for constructing this effect, for example, 
in the superfluid '^He Joesphson weak link \§\, high den- 
sity phases of QCD [9], high Tc superconductor system 
with SO (5) symmetry [10| and nonlinear optics induced 
non- Abehan field However, there are no any specific 
experimental constructions so far. 

How to design an experimental protocol to observe this 
novel effect in future experiments? To our knowledge, 



this effect has not been explicitly and simply spelled out 
in any real physical system, which is what we attempt 
to do here. In order to generalize to the non-Abelian 
junction in experiments, we need a system of multi- 
component order parameter which has a non-Abelian 
symmetry in the order parameter space. In contrast with 
magnetic tarp, the spin of the alkali atoms are essen- 
tially free in an optical trap [l^, [H, This spinor 
nature properly provide the scenario of our non-Ablian 
construction. In the following, we briefly introduce the 
system about a spinor atomic BEC in a double-well op- 
tical potential. Although the dynamical tunneling prop- 
erties of spin-1 and "pseudo spin-1/2" bonsonic systems 
were calculated [ll,|il,[l3,[li, 19], the essence of the non- 
Abelian effect has not been captured yet. In present Let- 
ter, we focus on the pseudo Goldstone modes due to the 
non-Abelian symmetry breaking, which is at the heart of 
Josephson effect. For concrete construction, we propose 
spin-2 BEC in double optical trap. The spin-2 system has 
possible advantages, compared to the spin-1 system, in 
the sense that the symmetry properties are much richer 
to explore non-Abelian effects. 

Ground state structure. — Let us consider a system of 
a homogenous spin-2 Bose gas with s-wave interaction. 
This system can be described by the following mean field 
free energy 
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J^W = ^[co(^W + ci(f)^ + ||ep] 



M^^^, (1) 



where tp = {ip2, tpi, tpo, fp-i, '0-2)"^ is the order parameter 
of the spin-2 Bose system, co,ci and C2 are interaction 
strengths related to the scattering length in different spin 
channels, (f) — ^p^itp is the mean value of the spin opera- 
tor and 9 = J2a=-2(~^)°"^a'^-a represents a single pair 
of identical spin-2 particles. The ground state configura- 
tion can be determined by minimizing this free energ; 
There are several distinct phases in this system [20, 2l| 
We will discuss these phases under zero magnetic field 
and analyze the symmetry and low-lying excitation spec- 
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trum of each phase. 

(I) Ferromagnetic phases: When ci < and ci — 
C2/2O < 0, two kinds of ferromagnetic phases are en- 
ergetically favored. The corresponding ground state 
configurations are given by -0 = ■\/ne'^(l, 0, 0, 0, 0) or 
-0 = V"e*^(0, 1,0,0, 0), where n = A*/(co + 4ci) is the 
particle density and 9 is an arbitrary global phase. It is 
obvious that these ground states have a, U{1) symmetry 
which leads to only one massless Goldstone mode. There- 
fore, two uncoupled system have a U{1) (g) U{1) symme- 
try. This symmetry will soft breaks into aU{l) diagonal 
symmetry when a weak coupling is applied. This pat- 
tern of symmetry breaking corresponds to an Abelian 
Josephson effect. The low-lying excitation spectrum of 
this state has been derived as ujk = \J (k[(-k + 2g4n) [20| , 
where = fc^ and (74 = co-|-4ci. We should note that this 
Goldstone mode will break into two modes when coupling 
is applied: one zero energy mode and one pseudo Gold- 
stone mode. This pseudo Goldstone mode has a finite 
but small gap and leads to a density mode fluctuation in 
d.c. Josephson current. 

(II) Antiferromagnetic phase: In the absence of mag- 
netic field, there is only one kind of antiferromagnetic 
phase when C2 < and c\ — 02/ > are satisfied. The 
corresponding ground state configuration is degenerate 
with respect to five continuous variables [2l[. Accord- 
ing to Goldstone's theorem, there always exist massless 
modes in the ordered state when a continuous global sym- 
metry is spontaneously broken. As a result, the above 
degeneracy of ground state will lead to five massless 
Goldstone mode. Some of these modes will lead to non- 
Abelian Josephson effect because the fiuctuation of dif- 
ferent components of the spinor BEG in the each trap are 
coupled in the corresponding equations of motion. Sim- 
ilar to the ferromagnetic case, we will show that each of 
these modes leads to a pseudo Goldstone mode. These 
pseudo Goldstone modes will give rise to five different 
kinds of Josephson currents with different frequencies. 

(III) Cyclic phase: When ci > and C2 > 0, the cyclic 
phase is energetically favored. The ground state config- 
uration is given hy ip = y/n{^-^,0, ^75"'^' ^t^^re 
n — fi/co is the particle density and the global phase 9±2 
and 6*0 satisfy + 9-2-29o = tt We find that there 
are two independent continuous phase variables in the 
cyclic ground state. As a result, there should be at least 
two Goldstone modes. We will show that these modes 
also lead to non- Abelian Josephson effect. 

Non-Ablian Josephson effects. — We will analyze the 
Josephson effect of a spin-2 BEG system in a double- well 
optical trap, as shown in Fig. 1. We assume that the 
energy barrier between the two wells is strong enough 
so that the coupling between the Bose gas in each well 
is very weak and the overlap of the ground state wave 
functions in left and right well (which we denote as </3i(x) 
and i^_r(x)) can be safely neglected. We will also use 
the single mode approximation which means we take the 
same mode function for all five spin components, this is a 
widely used approximation and it is valid when the spin 
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FIG. 1: The experimental schematic of a spin-2 Bose gas 
trapped in a double well with chemical potentials /^l of left 
and of right trap, which initially satisfy fiL ~ f^R ■ To 
drive Josephson effect, we add a small distortion Sjj. to hr. 



interaction is symmetric. Under these assumptions, the 
system can be described by the following potential 

Vcouple = Vi^L) + - Jii^lijR + I^JjVl), (2) 

where ipL and ipR are the order parameter of the 
Bose system in left and right well respectively, J — 

- J d^x(^J^(x)(-V2 + K;en(x))v5_R(x) is the coupling pa- 
rameter, and V{ip) — F{il)) has been defined in equa- 
tion (1) with a redefinition of the following param- 
eters: Ci is redefined as Ci J cPx\(pl{x.)\'^, ji is rede- 
fined as jj. = -/ d3x(/?2(x)(-V2 -I- K,eii(x))(/3L(x) = 

— / d'^x(/3jj(x)(— + Vweii{y^))VR{^) which is just the 
single particle energy in each well. It should be noted 
that we have taken the same chemical potential for the 
Bose gas in right and left well, because we will only be 
interested in the d.c. Josephson effect which captures the 
essence of the non- Abelian symmetry breaking as simple 
as possible. With this potential, we can derive the equa- 
tion of motion of the Josephson current in each phase and 
analyze the relationship between the coupling parameter 
and the pseudo goldstone modes. However, as we have 
mentioned above, the symmetry of the ground state in 
ferromagnetic phase only leads to an Ablian Josephson 
effect which is not interested in present paper. There- 
fore, we will just analyze the antiferromagnetic phase 
and cyclic phase which are important realizations of non- 
Ablian Josephson effect. 

(I) Antiferromagnetic phase: To obtain the equa- 
tion of motion of the Josephson current, we can add a 
small fluctuation around the ground state "0^ = 0^ = 
^(1, 0, 0, 0, 1) and we should first determine the ground 
state particle density n [l^ . By minimizing the potential 
Vcoupie under the above antiferromagnetic ground state, 
we get n = {^ + J)/{co + C2/5). With this ground state 
configuration, we can obtain the equations of the fluctu- 
ations in this phase and analyze the excitation spectrum. 

The m = mode. The equation of motion of m = 
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mode is given as 



d C2 C2 

«^0OL = + J)4'0L + Y^-'/'Sl - J4>0R, (3) 

and a similar set of equations of (t>oR- Since this mode 
is decoupled from others, it corresponds to an Abelian 
Josephson current. By solving equation ([3]), we ob- 
tain the eigenenergies of this mode, one is zero corre- 
sponding to a massless Goldstone mode and the other is 

ujo = 2^ J{J + ^l^^al" ^ corresponding to a pseudo Gold- 
stone mode. 

The m = ±1 coupled mode. The m — ±1 mode are 
coupled in the following equations 



[n{ci - ^) + J^iL 
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n{ci - —)<i)-iL - 
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dt 
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+ n{ci - —)(t)*iL - J4>-lRi 
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(4) 



and a similar set of equations of <f>±iR. The solution in- 
volves one zero energy oscillation and one pseudo Gold- 
stone mode with a gap oi u}±i = 2^n{c\ ~ ^) J -I- J^. 

The m = ±2 coupled mode. The m = ±2 mode are 
coupled in the following equations 



{ — + J)<P2L + — (P2L - 

1 2 

-7l(co - 4ci + -C2)(0-2L+C2l). 



{ — + J)<P-2L + — (, 



-2L 



2R 



1 2 

+ T^l'^ico - 4:Ci + -C2){(j>2L + <p2L), 

2 5 



(5) 



and a similar set of equations of <j)±2R- Through this 
set of equations, we obtain two zero energy modes and 

two pseudo Goldstone modes with energy gap uj^^\ = 

2^n{co + f)J + J'' and J^l = 2^n{ci - + .P . 
In Fig. 2 we give the dependence of the above frequencies 
on coupling parameter J. Recently, a polar behavior has 
been observed in the F=2 ground state of ^"^Rb conden- 
sate [1^ . We expect that the above modes of fluctuations 
can be observed in this system in future experiments. In 
the case of system, the value of interacting strengths 
under typical experimental condition are given as |22| : 
cin : — lOnK, C2n : — 0.2nK and cqti about 150nK. 
According to the weak coupling limit, we assume that the 
coupling parameter J is much smaller than the interac- 
tion energy of the condensate and given as about O.lnK. 
Under these conditions, we can obtain the frequencies of 
the fluctuation related to the antiferromagnetic phase, 
which is of order lOOHz. The measurement of fluctua- 
tions on this characteristic time scale (about lOnis) is 
accessible in current experiments. 
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FIG. 2: The frequencies of the pseudo Goldstone modes as a 
function of coupling parameter J in the case of antiferromeg- 
netic phase. All the frequencies is proportional to when J 
approaches zero and the dependence on J becomes linearizing 
when J is large compared to the interaction energy. 



(II) Cyclic phase: Following the same procedure in the 
antiferromagnetic phase we can get the equations of mo- 
tion for the cyclic phase. 

The m = ±2,0 coupled mode. The fluctuations of 
m — ±2 and m = mode are coupled together and can 
be described by the following equations 



( ^ + J)02L + ^— 02L 

a — f3 1 a ,* 

<?^-2L + 4 (a- /3+7)'/'-2L + ^-/-OL 

a — 7 



2V2 

{^ + J)c^,L + ^{<t>2L+<t>^2L) 

■^4>Il + :^{4'*2l+4'*^2l) - J(t>OR, 

Q! + /3 + 7 ^ + f^A* 



+ ^02L + J(a-/3 + 7)02L + ^</'5L 



a — 7 
'2V2^ 



lOL ~ J4>-2R, 



(6) 



and a similar set of equations of (j)±2R and 4>qb.- The 
parameters in the above equations are defined following, 
a = cqu, l3 — Acin and 7 = 2c2n/5. By solving the 
above equations, we find that each Goldstone mode in 
the corresponding uncoupled system (20j break into one 
massless mode and one pseudo Goldstone mode. Since 
there are two Goldstone modes in the uncoupled system. 
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we find two pseudo Goldstone modes with energy ojq = 

2VJ2 + aJ and Jf^l^ = 2^ + f J. 

The TO = ±1 coupled mode. The fluctuations of to = 
±1 mode are also coupled together and can be described 
by the following equations 

^-J^9lL = (^ + J)(PlL + ^<P-lL + ^^<PlL~J<PlR, 
«^'/'-lL = (-2 +'^)'P-li+4'?'lL ^<P-lL--''/'-lfl.(") 

and a similar set of equations of (j)±iji. As we know, there 
are two massless Goldstone modes with the same energy 
in the uncoupled system ^2Q] . By solving the above equa- 
tion of motion, we can see that each of them leads to a 

pseudo Goldstone mode with a gap uj±i = 2yJ ,P + ^ J. 
These kind of fluctuations in cyclic phase is expected to 
be realized in a condensate of ^^Rh atoms [21J. Based on 
the current estimates of scattering lengths, the value of 
interacting strengths under typical experimental condi- 
tion are given as: Cin : — 20nK, C2n : — 0.6nK and 
cqtl about 600nK. Under these conditions, we can also 
estimate the d.c. Josephson frequencies in cyclic phase 
which is about 100-300Hz. The dependence of the above 
frequencies on coupling parameter J is shown in Fig. 3. 

All the above analysis is base on the assumption that 
the phase of the condensates in two wells are equal to each 
other in their ground state. However, there is another 
experimental achievable configuration called tt state in 
which the phases of condensates in right and left well 
have a tt difference. We find that, analogue to spin-1 and 
single component BEG system [T^, [H, HSIj the oscilla- 
tion around this state could be unstable in some cases. 
For example, in the case of Antiferromagnetic phase, 
the value of pseudo Goldstone mode t^o will change to 

t^o - 2y^J(J- ^Ifi^ in TT state. When J < this 
mode can clearly become unstable. This in fact happens 
in all the pseudo Goldstone modes we calculated above. 
As a result, the experimental protocol we propose below 
should only restricted in the oscillation around the true 
ground state of the system not around the tt state. 

Experimental signatures of Non-Abelian Joesphson 
effects. — We want to propose an experimental proto- 
col to observe the non-Abelian Josephson effect. The 
experimental set up of a spin-2 Bose gas trapped in a 
double well is illustrated in Fig. 1. The d.c. non-Abelian 
Josephson current can be detected with the following 
steps. The first step is to initiate a density oscillation 
in the system. This can be realized by slightly changing 
the depth of one well, which will cause a small imbalance 
in chemical potential {fXR + 6fi) between the two 

wells, and then tune it back. The next step is to detect 
the time dependence of the particle numbers in differ- 
ent spin component. Such kind of detection can be real- 
ized by first spatially separating different spin component 
with a Stern-Gerlach method during time of flight after 
switching off the trapping potential. Then, the number 
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FIG. 3: The frequencies of pseudo Goldstone modes as a func- 
tion of coupling parameter J in the case of cyclic phase. The 
dependence on J is similar to the case in antiferromegnetic 

(2) 

phase. However, there are two degenerate modes uJq j_2 and 
uj±i in this case, which are originated from the two degenerate 
Bogoliubov modes in the uncoupled system. 



of atoms in each spin component will be related to the 
respective spatial density distributions which can be eval- 
uated by the absorption imaging method. Following the 
above steps, one can measure the density oscillation in 
each spin component which are coupled together due to 
the non-Abelian symmetry of the system. The measure- 
ment in the dependence of the oscillation frequencies on 
J can be realized by varying the barrier between the two 
wells and repeating above measurement. 

In mean field theory, the condensates of atom in 
F=2 state is predicted to be polar (ci — C2/2O > and 
C2 < 0), but close to the border to the so-called cychc 
phase (ci > and C2 > 0) [2l|. Furthermore, polar be- 
havior in the F=2 ground state of ^"^ Rb has been observed 
in recent experiment [2^ . As a result, we expect that the 
pseudo Goldstone modes of the antiferromegnetic phase 
could be observed in experiments. As we have men- 
tioned, the value of interacting strengths under typical 
experimental condition are given as [22l |: cin : — lOnK, 
C2n : — 0.2nK and cqu about 150nK, which leads to 
the time scale of about 10ms in the fluctuations in this 
system. On this time scale, the measurement we pro- 
posed above is completely accessible in recent experiment 
in F=2 spinor Bose-Einstein condensates of system 
[2^ . Although there is still no such kind of measurement 
performed in a system with cyclic phase, we expect that 
it will be realized in a condensate of ^^Rb atoms in near 
future [21]. 

In summary, we propose an experimental protocol to 
realize the so called non-Abelian Josephson effect in a 
spinor Bose system. We find that the frequencies of 
pseudo Goldstone modes do not only relate to the cou- 
pling parameter but also to the interacting strengthes, 
which is a nonlinear effect due to the spin dependent in- 
teraction. Our results are of particular significance for 
exploring the new features of the non-Abelian Josephson 
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